X-ray absorption spectra at the Ca-L2,3-edge calculated within multi-channel multiple 

scattering theory 

Peter KriigeiQ 

LRRS, UMR 5613 Umversite de Bourgogne - CNRS, B.P. ^7870, 21078 Dijon, France 



o 
o 



(N 



C/2 



Calogero R. Natoli 
INFN Laboratori Nazionali di Frascati, Casella Postale 13, 1-00044 Frascati, Italy 

(Dated: February 2, 2008) 

We report a new theoretical method for X-ray absorption spectroscopy (XAS) in condensed matter 
which is based on the multi-channel multiple scattering theory of Natoli et al. and the eigen- 
channel R-matrix method. While the highly flexible real-space multiple scattering (RSMS) method 
guarantees a precise description of the single-electron part of the problem, multiplet-like electron 
correlation effects between the photo-electron and localized electrons can be taken account for in a 
configuration interaction scheme. For the case where correlation effects are limited to the absorber 
atom, a technique for the solution of the equations is devised, which requires only little more 
computation time than the normal RSMS method for XAS. The new method is described and an 
application to XAS at the Ca L2,3-edge in bulk Ca, CaO and CaF2 is presented. 

PACS numbers: 71.15.Qe, 78.70.Dm 



X3 

o 



> 

00 
(N 

O 

o 



I 

o 
o 



X 



I. INTRODUCTION 

Multiple scattering (MS) theory provides an accurate 
and flexible scheme for the calculation of unoccupied elec- 
tronic states which are probed by various synchrotron ex- 
periments such as X-ray absorption spectroscopy and res- 
onant elastic and in-elastic X-ray scattering. The stan- 
dard theory relies on the single-particle picture, that is, 
it neglects electron correlation effects. This is a great 
shortcoming, since core-level X-ray spectra are often 
strongly modified by electron correlation, in particular 
by the Coulomb and exchange interaction of the valence 
electrons with the core hole. In transition metal and 
rare earth systems, this interaction can give rise to pro- 
nounced atomic multiplet and satellite structures in the 
spectra, which can only be accounted for through many- 
electron calculations. A generalization of MS theory to 
many-electron wave functions was developed by Natoli 
et al. [1] and is known as "multi-channel" MS theory. 
Probably the most difficult part of this approach is the 
calculation of the intcr-channcl potential. Here, we pro- 
pose a reformulation of the theory, where the latter prob- 
lem is completely avoided. Instead, the multi-channel T- 
matrix is calculated variationally using the eigen-channcl 
R-matrix method H,|1,^|EI1|' While R-matrix methods 
are well known in atomic spectroscopy, they have, to our 
knowledge, never been used for condensed matter prob- 
lems. Michiels et al. 0] presented a calculation of electron 
energy loss from NiO using an R-matrix method. They 
used, however, an atomic model where all solid state ef- 
fects were described phenomenologically using an crystal 
field and a reduced Coulomb interaction. 

Here, we present a new formalism for X-ray absorption 



in condensed matter, based on multi-channel MS theory 
and the eigen-channel R-matrix method. It allows to take 
account for local electron correlation effects in a multi- 
channel, that is, configuration interaction scheme. At 
present, the type of correlations that can be handled on 
this level are limited to those between one electron in a 
delocalized state and a finite number of electrons/holes 
in (sufficiently) localized orbitals. Sufficiently localized 
means that the wave function is negligible small beyond 
the atomic radius. This applies exactly to inner-core 
shells and well to the 4/-shell in rare earths. Extensions 
of the method to include correlation effects between sev- 
eral delocalized electrons are under way. 

In this paper, we present the formalism and report 
results on the Ca-L2,3-edge absorption of diflFerent Ca 
compounds. The Ca-L2,3-edge is an interesting test case 
for the new method, because the L2 and L3 absorp- 
tion channels are strongly coupled through the photo- 
electron - core-hole Coulomb interaction. This leads to 
a branching ratio of about 1:1, far from the statistical ra- 
tio (2:1) which is obtained in single-particle theory. From 
a point of view of atomic multiplet theory |^ 0, |^ , the 
non-statistical branching ratio is easily understood as a 
case of strong intermediate coupling in the {2p^3d^) fi- 
nal state. The multipole and exchange part of the 2p- 
3d Coulomb interaction (Slater-integrals F'^, with 
/c > 0) is of comparable strength as the 2p-spin-orbit 
interaction, which gives rise to correlated 2p^id^ final 
state wave functions, where the 2pii2 and 2p3/2-holes 
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are strongly mixed. The fact that the branching ratio 
does not change when going from atomic Ca [sl to var- 
ious Ca compounds 11] is empirical evidence that this 
atomic multiplet picture remains valid in condensed mat- 
ter. However, a purely atomic model is not sufficient to 
account for fine structure in the L2,3-edge spectra, which 
depends strongly on the atomic environment pT| (and 
which is thereby of practical importance for structural 
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and electronic analysis.) Atomic models including crys- 
tal field have proved quite successful in rep roducing the 
experimental spectra at the L2,3-edge In that 

approach, all extra-atomic effects are, however, treated in 
an empirical way, by introducing adjustable parameters 
for crystal field and (possibly) band broadening. Zaanen 
et al. was the first to go beyond the atomic model by 
considering a model Hamiltonian that included not only 
the atomic 2p-3d multiplet coupling but also the single 
electron density of states of bulk Ca. The electron-hole 
problem was solved exactly using a Green's function tech- 
nique. While being physically sound, Zaanen's method 
was not fully based on first principle calculations, but 
introduced a number of empirical parameters. Later on, 
Schwitalla and Ebert calculated the spectra in the 
time-dependent local-density-approximation (TD-LDA) . 
For bulk Ca, they obtained the correct branching ratio, 
but the fine structure of their spectrum was quite differ- 
ent from the experimental one. Recently, Ankudinov et 
al. 10 studied the branching ratio problem with a gener- 
alization of TD-LDA. By adding a frequency and matrix- 
element dependent exchange-correlation contribution to 
the TD-LDA kernel, they obtained a branching ratio in 
good agreement with experiment for Ca and the whole 
transition metal series, while in Ref. |0| this was true 
only for the lighter elements (from Ca to V). From the 
theoretical studies cited above, it may seem that the 
branching ratio problem at the i2,3-edge of Ca has been 
thoroughly investigated. Despite of this, we have chosen 
the Ca system as a test case for our new method, which, 
we believe, provides new insight into other aspects of the 
problem, like the orbital relaxation around the core-hole 
and the reason for the need of a 20% reduction of the 
Slater integrals F'^ and G'^' in atomic multiplet calcula- 
tions Jl0| . The novelty of the present method comes in at 
two levels: it is the first application of the multi-channel 
MS formalism and it is (to our knowledge) the first true 
application of i?-matrix techniques to a condensed mat- 
ter problem. The combination of these two features will 
allow us to shed some light on the two points mentioned 
above (orbital relaxation and reduction factors) and to 
present an application to the Ca compounds CaO and 
CaF2, in which ligand field effects and multiplet struc- 
ture are treated in a unique framework in an ab initio 
way. 



The paper is organized as follows. In section ^ the 
more general aspects of the formalism are outlined. Fur- 
ther details about the multi-channel MS theory can be 
found in the appendix. In section IIIII the formalism is 
applied to the L2,3-edge absorption of 3d" systems with 
an emphasis on the screened electron-hole interaction in 
the final state. In section IWI some numerical aspects are 
discussed. In section Ivl results are presented for bulk Ca, 
CaO and CaF2. Finally, some conclusions are drawn in 
section IVTl 



II. GENERAL FORMALISM 

In the present approach we go beyond the independent 
particle model by considering a correlated wave function 
for a finite number of N electrons. All other electrons are 
described within the independent particle approximation. 
Among the N explicitly treated electrons, at most one 
is in a delocalized orbital, all others necessarily occupy 
localized orbitals. By definition, a localized orbital is 
one that is negligibly small outside the atomic sphere. 
This applies exactly to inner shell orbitals but also to a 
good approximation to 4/-orbitals of the rare-earths. In 
the ground state wave function, the N electrons include 
the core electron that is excited in the XAS process plus 

— 1 other electrons in localized orbitals. The XAS 
final state wave function then contains A'^ — 1 localized 
electrons and one electron (the "photo-electron") in a 
delocalized state above the Fermi level. In other words, 
we consider a correlated final state wave function that 
couples the photo-electron with the core- hole and and/or 
a finite number of other localized electrons. 

In order to make the derivation less abstract, we shall 
now consider the specific case of XAS at the Ca-L2,3 
edges. The formulae are kept general and can easily 
be applied to other systems to be described with corre- 
lated wave functions satisfying the above requirements. 
For the ground state we consider the six electron wave 
function made of the 2p core electrons. The initial state 
with energy Eg is thus simply given by the closed 
shell configuration {2p^,-^So). Final states have energy 
E — Eg + hu! and a (2p^e^) configuration, where e de- 
notes a (one-electron) state in the continuum above the 
Fermi energy. The crucial point is that we take into ac- 
count multiplet effects through a configuration interac- 
tion ansatz for the final state wave function, which is 
developed as 

^ = AY,^a{X)cf>^{x) . (1) 

a 

Here <&q. is one of the six (2p^) states, labeled by a— 
Ucfic) {jc = 1/2,3/2, = ~jc---jc); X collects ah 
core-electron coordinates. The (2p^) multiplet energies 
are Ea = Eg — edjc), where edjc) are the negative bind- 
ing energies of spin-orbit split 2p{jc) levels. For each 
$a , there is a component (j)a of the photo-electron wave 
function. The (radial, angular, and spin) coordinate of 
the photo-electron is denoted x = (r, x,cr). Finally, A 
denotes the anti-symmetrization operator. 

Multi-channel multiple scattering. The total photo- 
absorption cross section is calculated using the multi- 
channel multiple scattering method by Natoli et al. P|. 
As shown in detail in the appendix, it is given by 

ct(w) cx wim 1^ M^T^^,Mr' | . (2) 

Here r=aLs is the set of all quantum numbers of 
with L = Im being the orbital and s the spin quantum 
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numbers of the photo-electron. Mr = (^'™|D|^'g) are 
the transition matrix elements; we consider only dipole 
transition in the length approximation. is the inside 
solution that matches smoothly onto the outside solution 

vp-t ^ J2 <f r' iXxa)Zr'r{r)/r . (3) 
r' 

Here, we have introduced <i>r = 'ta{X)YL{x)Ss^a-- "In- 
side" and "outside" refer to the atomic sphere of the ab- 
sorber, i.e. r < vq and r > tq, respectively, vq being the 
muffin-tin radius. The matrix Z of radial photo-electron 
functions is given by 

Zrr'{r)/r = ji(fcQ^)[^o ^]rr' - iKhf {kar)5YT' ■ (4) 

Here, = ji + ini and ji,ni are the usual spherical 
Bessel and Neumann functions, fc^ is the wave number 
of the photo-electron, given by k'^ + Vq — €a = E — Ea, 
where Vq is the interstitial potential, ^rr' is the multi- 
channel atomic T-matrix of the absorber (at site i—0). In 
Eq. ^ , Tpp, is the multi-channel scattering path operator 
connecting sites i and j. It is calculated for a finite cluster 
by inversion of the matrix m = T~^, whose elements are 
given by 

TOrr' = Sij[t~\r' - Saa'kaG''l,^,{ka)Sss' . (5) 

Here, ti is the multi-channel atomic scattering matrix 
of atom i, and G^^, are the real space KKR structure 
factors. 14] Apart from the absorber, we treat all atoms 
in the standard onc-clcctron muffin-tin approximation, 
which implies ti,rr' — ^i;(^a)<5rr', for all i^O. Since 
these T-matrices for iy^O as well as the structure factors 
G^Ij^, are single-channel quantities, the only channel-off- 
diagonal terms of m are located in i=j=0 block. This 
particular structure of the m-matrix allows us to use an 
efficient partitioning technique for the inversion of m. 

Partitioning technique. We divide the system into 
absorber atom {i=0) and "environment", i.e. all other 
atoms with i^O, collectively labeled 'e'. For the absorp- 
tion cross section we need only the absorber block of 
the r-matrix. Using simple matrix algebra, this quantity 
can be expressed as 

In the second equality, we have used m^^ — ^o^ 
introduced the reflectivity p = mP'^[Trf'^]^^irf^ , which 
contains all the information we need from the environ- 
ment. Once p is known, the remaining problem is a 
purely atomic one. Now p is diagonal in the channel 
indices a since it does not involve the i=j=0 block of the 
TO-matrix. It can therefore be calculated using standard 
(single-channel) MS theory. Explicitly, we have 

Prr' = Saa'PLL'{ka)6ss' , (7) 

where 

PLL'ik) = k'J2(^L\4k)r%^,,,{k)Gi",,,^,{k) . (8) 



Here, the sum runs over L", L'" , ij^O, j^O, and f is the 
single-channel r-matrix of the system without absorber 
(a-diagonal terms of [m'^'^]~^). 

Eigen- channel R-matrix method. The remaining prob- 
lem is the calculation of the multi-channel T-matrix of 
the absorber and the inner solutions ^f™. This is done us- 
ing the eigen-channel i?- matrix method |^ ^ |^ . In 
the following we recall some basic features of this method 
for the convenience of the reader and in order to intro- 
duce our notation (which follows most closely that of 
Ref. 0). The R-matrix is a multi-channel generaliza- 
tion of the logarithmic derivative of the radial wave func- 
tion. As reaction volume, we use the atomic (or "muffin- 
tin") sphere of the absorbing atom with radius tq. With 
Eq. (jSJ, the R-matrix can be defined as 

y^^Rrr"Zr"r'{ro) = Zrr'jro) . (9) 
r" 

Here we have introduced the notation X = dX/dr. Using 
Eq. Q and its derivative with respect to r, the t-matrix 
can be readily calculated from the R-matrix as 

t-^ ::=iK{Rj - J)-\RH - H) . (10) 

Here all the quantities are matrices with indices TV and 
are evaluated at r — tq. Furthermore, the quantities 
K, J, H are diagonal matrices with elements K^r — k^, 
Jrr = karojiikaro), and Hrr = kar^hl {karo). 

In the eigen-channel method, the R-matrix is obtained 
directly in diagonal form; for given energy E, a basis of 
eigenstates and eigenvalues bk is found^ysolving the 
following generalized eigenvalue problem j5|, |6j . 

{E-H-L)^k = Q^khk. (11) 

Here H is the Hamiltonian, L = X^iLi ^{'''i ^ fa)-^-^Ti 
is the Bloch operator that restores Hermiticity of H in 
the finite reaction volume, i.e. the atomic sphere and 
Q = Y^f^i5{ri — Tq) projects onto its surface. Among 
all solutions of Eq. 1)11(1 . only those with \bk\ < oo are 
physically acceptable. Their number equals the number 
of channels F In order to solve Eq. Hll|l we develop 

*fe = XI "^rucru^k (12) 

with trial functions of the form 

*r. = A{'i>r{Xxa)P,ir)/r} . (13) 

As radial basis functions P^, we use solutions of the 
radial Schrodinger equation for angular momentum I 
and a spherically symmetric, local one-electron poten- 
tial V(.f{(r). In the present application, we take for Vcff 
the sum of the ground state potential Vg and a partially 
screened core-hole potential Vc (see Eq. (|16|) below). As 
usual in the eigen-channel method, we use closed-type 
orbitals with boundary conditions Pi^iro) = 0, and open- 
type orbitals with boundary conditions dP^/dr{rQ) = 
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0. Since 2p es transitions have neghgible intensity 
in the near-edge region, we here include only I = 2, 

1. e. c?-waves in the basis. The generalized eigenvalue 
problem in Eg. Ullfl is solved using standard numeri- 
cal routines The eigenvectors of the R- matrix are 
given by Wrk = VNi^q J $r^'fc, where the integration 
is over Xxa and the remaining radial coordinate of '^k 
is taken at 0. The factor \/N comes from anti- 
symmetrization. From the orthogonality of the channel 
functions <i>r and Eqs H12I13() we have 

Wrk = ^cri/,fcP,.(ro) . 

We normalize the generalized eigenvectors cri,,k {k fixed) 
such that J2r l^^r/cP — 1- Then W is unitarian and the 
R-matrix is given by 

Rrv ^ -J2^Tkb^'wl^, . 

k 

The inner solutions that match the outer ones are 

E *rvCrv.feW^2r"^r"r(ro) . 
r'r"uk 



in Eq. which we denote as iTiy) = \2p^vd^,T) with 
T—jciJic'fns. We have 

{Tv\Ha\T'v') = [Eg - edjc) + e„)Srr'S^„' , 

where 

S,,, EE / drP,{r)P,,{r) 
Jo 

is the overlap integral, is the energy of the P^, orbital, 
and the other quantities have been defined before. Note 
that 6rr' is ensured by the orthogonality of the angular 
and spin functions. For the calculation of the matrix el- 
ements of V, we make a basis transformation from the 
uncoupled states \2p^jriJ,r., vd}ms) to LS coupled states 
\2p='v(f, {LS)JM) lig. In the LS coupled basis, the ma- 
trix elements of V are given by 

(Fi^lFlrV) = [wC'^+'QUdrr' , (14) 

where now T={LS)JM. The w's can be expressed in 
terms of the following generalized Slater integrals 

/••To pro 2j-k 

Fl,^ dr dr'P^^{r)PAr')^P2p{r)Pu'{r'). 
Jo Jo ?'> 



III. ELECTRON-HOLE INTERACTION 

We describe the sub-system of N electrons through a 
Hamiltonian of the form 

N 

H^^'> =H„ + V = Y,ho{^) + V . 

i=i 

Here ho is the one particle Hamiltonian of the chosen in- 
dependent electron model and i is an electron label. If N 
was the number of all electrons in the system (iVaii), the 
exact perturbation V would be given by the bare two- 
particle electron-electron interaction terms minus the ef- 
fective electron-electron potential VeS that is included 
in ho- However, since in our case TV ^ A'aii, there is 
no (simple) exact expression of iJ^^^ and the "best" ap- 
proximation for V is not necessarily given by the exact 
expression of the case N — iVaii. The reason is that VcB 
and thus ho are determined by iVaii rather than only N 
electrons, and the Coulomb interaction in V is screened 
by the A^aii — N other electrons. 

For the system studied here, these considerations are 
of interest only for the final state. The ground state, be- 
ing a closed shell configuration {2p^,^S), is well described 
by a single Slater determinant with the {2p) orbitals cal- 
culated from ho- For the {2p^e^) final states, the per- 
turbation V is the screened photo-electron - core-hole 
Coulomb interaction. 

We shall first take for V the unscreened interaction and 
discuss the effect of screening below. We have to calculate 
the matrix elements of H, L and Q for the basis states 



fro fro 

Gt,= dr dr'P2p{r)P,{r')^P,,{r)P2p{r') . 
Jo Jo 

Here, f>(<) is the larger (smaller) of r and r' . The ex- 
pressions for the wi^^^^L) are given in Ref. [T7|: 

w{^P) = -F" - F'^/5 + 40749 

w{^p) = -F" - pyb 

w{^'^D) = -F° + Fy5 
w{^F) = -F^ - 2F^/35 + I8G749 
w{^F) = -F" - 2F^/35 
where, in our case, all quantities are matrices with in- 
dices {vi'')- The matrix elements of the other operators 
needed in the eigen-channel method are easily calculated: 

{Tv\E\T'v') ^EStt'Su^' , 
{Vv\Q\V'u') ^ 5TT'PAro)P.'{ro) , 

(r^.|L|rV') =<Jrr'P.(ro)^(ro) . 

The dipole operator selects (^P) basis states and thus 
only J — 1 final states give a contribution. The reduced 
matrix elements are non-zero for r=(^'^+^Li,A/=l) and 
given by 

(vE'.llrllM'i?) = -2 ^ hpJroMr'Zr'Tiro) (15) 

vkF' 

with ro=(iPi,Af=l) and hp,^ = J^" P2p{r)rPAr)dr. 
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Screening model. As it is well known from cluster and 
impurity model calculations the monopole term of 
the electron-hole Coulomb interaction (corresponding to 
the Slater integral F") is drastically screened, while the 
higher order multipole and all exchange terms (Slater in- 
tegrals , with fc > 0) are essentially unscreened. 
Let us note that in multiplet calculations also the higher 
order terms are generally reduced from the calculated 
values. 1^ Il0| The need for this reduction of some 20% is, 
however, not due to screening, but comes mainly from 
the neglect of configuration interaction in the single- 
configuration multiplet approach As will become 
apparent in the next section, the relevant configuration 
interaction is included in our approach, so that there is 
no need for reduction of the Slater integrals with fc > 0. 

We therefore apply screening only to the monopole 
term 2/r> of the Coulomb operator 2/|x — x'|. This 
defines the (unscreened) multipole part V^ = 2/|x — x'| — 
2/r> of the interaction V . In the space of trial func- 
tions we have chosen, the operator 2/r> is diagonal 
in r. Within this space, it is therefore equivalent to a 
one-electron potential Vu{t), namely the Hartree poten- 
tial of a spherically symmetric core-hole which is given 
by : Vu{r) = J dr'[P2p{r')Y /r-^. We can thus handle 
screening of the monopole term on a single-particle level 
by replacing the unscreened core- hole potential Vu{r) by 
a screened one Vc{t), which we add to /iq- In this way the 
effective potential VeS used in will not be the ground 
state self-consistent potential Vg but VcS = Vg + Vc. 

A simple approximation for Uoff is given by the fully 
statically screened potential Wsuporccii which is obtained 
from a self-consistent super-cell calculation with a core- 
hole on the absorber site. This effective potential, which 
features full orbital relaxation around the core-hole, is 
frequently used in single-electron XAS calculations. We 
shall denote the corresponding core-hole contribution by 
Vs, i.e. Vs = i^suporceii — Vg. As will bccomc apparent in 
the result section below, the line shapes obtained with Vg 
are not satisfactory. We shall therefore allow for incom- 
plete screening by using a linear mixture between the 
unscreened core-hole potential Vu and the fully screened 
one Vs : 

Vc{r) ^ avu{r) + {1 - a)vs{r) , (16) 

where a G [0, 1] is an empirical parameter. As can be 
seen from the results below, a value of a « 0.1 gives 
best agreement with experiment. This fact indicates that 
orbital relaxation around the core-hole is overestimated 
in UsupcrccU which, we recall, is obtained from a super- 
cell calculation in the local density approximation (LDA) . 
This finding was to be expected, since it is known that 
in LDA the self-interaction of an electron is not exactly 
compensated as in the Hartree-Fock scheme, giving rise 
to over-relaxation and band gaps that are systematically 
too small compared to experiment. Probably the same 
calculation with Self-Interaction Corrections would cure 
this drawback. This will be the subject of a future in- 
vestigation. In the meantime we regard a as useful pa- 



rameter describing the correct amount of relaxation. Let 
us also note that an a priori estimate of this quantity 
could be obtained in a multi-channel MS theory that 
starts from the fully relaxed state and mixes in very many 
charge transfer excitations. This possibility shall also be 
explored in the future. 

In summary, the present treatment of screening con- 
sists in (i) replacing the screened electron-hole Coulomb 
interaction V by its unscreened multipole part V — 
2/|x — x'l — 2/r> and (ii) adding the partially screened 
core- hole potential Vc{r) in Eq. H16|) to the single-particle 
Hamiltonian ho. Point (ii) results in a modification of 
all radial wave functions Pv{t) and corresponding ener- 
gies e^, whereas point (i) simply removes all monopole 
terms (F^^,) from the interaction matrix H14|l . 



IV. NUMERICAL ASPECTS 

The standard MS calculation for the reflectivity of 
the environment Prj/jk) has been performed using the 
CONTINUUM code [l^. Finite clusters containing at 
least nine nearest neighbor shells around the absorber 
were used for all systems, such that the XAS spec- 
tra were well converged with respect to cluster size. 
The effective single-particle potential was calculated self- 
consistently in the local density approximation using the 
linear- muffin-tin-orbital method ^0]. In all systems we 
used space-filling (and thus partially overlapping) atomic 
spheres. In the compounds CaO and CaF2, we chose the 
relative atomic radii in such a way that the potential 
value on the sphere was approximately equal for Ca and 
the ligand, while keeping the overlap volume small. For 
CaO, a ratio of 3:2 between the Ca and O radii was found 
appropriate. For CaF2, the insertion of one empty sphere 
(E) per formula unit was necessary to keep the overlap 
small. We chose the ratio of the sphere radii of Ca:0:E to 
be 6:4:5, approximately. The fully screened core-hole po- 
tential Vs was obtained from super-cell calculations with 
a (spherically symmetric) 2p-hole on the absorber atom. 
We used a 32 atom simple cubic super-cell for Ca metal, 
and an fee 2 X 2 X 2 super-cell for CaO and CaF2. We 
found that the core-hole has only a small effect on the po- 
tentials of the neighboring atoms and that, consequently, 
it makes hardly a difference for the spectra whether the 
reflectivity is calculated with or without the core-hole po- 
tential. On the absorber atom, however, Vc{r) is strong 
and has a dramatic effect on the line shape as will become 
apparent below. 

For reasons of numerical stability, the reflectivity was 
calculated at complex energies with a small imaginary 
part, such that the spectra are effectively broadened with 
a Lorentzian function of about 0.3 eV FWHM. In order 
to simulate finite experimental resolution, the spectra in 
the result section were further broadened with a Gaussian 
function of 0.3 eV FWHM. 

In the eigen-channel method, convergence with respect 
to the number of radial basis functions has to be achieved. 
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il FIG. 2: CaO spectra obtained with the minimal basis set 

S; V/j ii'j (nc, no) = (1,1), along with the converged one (c). The latter 

f y y \ I was multiplied by a factor 1.4 for easy comparison. Spectrum 

V (b) was calculated with the Slater integrals F'',G'' reduced 

. . . . . . . . I by 20%. 
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FIG. 1: Convergence of the spectra with respect to the num- 
ber of radial basis functions in the example of CaO. The num- 
bers of closed-type (nc) and open-type (no) basis functions are 
indicated as {nc,no). 



In Fig.^we show this convergence in the example of CaO. 
The different basis sets are indicated as (nc,no), where nc 
(no) is the number of closed-type (open-type) functions. 
We start from functions without nodes in < r < rg and 
increase the number of nodes one by one. For example, 
the (3,1) spectrum was obtained with three closed type 
functions of zero, one and two nodes and one open-type 
function of zero nodes. Figure ^ shows the converged 
spectrum with basis set (6,2) in the upper panel and dif- 
ference spectra with respect to (6,2) in the lower panel. 
It can be seen that five closed-type and only one open- 
type function are sufficient for good convergence. For 
the spectra in the results section below, we have used 
the (5,1) basis set. 

It is interesting to note that one can considerably im- 
prove the spectrum calculated with the minimal (1,1) ba- 
sis set by reducing the values of the Slater integrals F'', 
G'^ (artificially) by some 20%. Figure [21 shows the (1,1) 
spectrum with full (a) and 20% reduced (b) values of 
the Slater integrals, along with the converged spectrum 
(c), which was multiplied by 1.4 for easy comparison of 
the peak ratios. Clearly, a 20% reduction of Slater in- 
tegrals improves considerably the (1,1) line shape, both 
as far as peak positions and relative peak intensities are 
concerned. Apart from the overall amplitude, which is 
about 40% too big, the spectrum almost coincides with 
the converged one. This result is closely related to the 
fact that in atomic single configuration multiplet calcula- 
tions, reduction factors of 10-25% for the Slater integrals 



are generally needed to make the relative multiplet en- 
ergies and line strengths agree with experiment |lOl Il8| . 
Such a rescaling procedure effectively accounts for con- 
figuration interaction that lies beyond the single config- 
uration calculation 18], namely coupling to higher lying 
electronic configurations. Precisely this feature is seen in 
Fig. 13 when one realizes that in the minimal set (1,1) de- 
scribes essentially only the 3d orbital, while in the (5,1) 
basis set of the converged spectrum, all nd orbitals up to 
n = 7 are included. 

In the practical implementation of the method, we first 
calculate the reflectivity matrix pLL'{k) on a fine mesh in 
the relevant (photo-electron) energy interval. In a second 
step the atomic multichannel calculation is performed 
for each total energy E = Eg + uj. The R-matrix and 
the inner solutions ^'P are calculated through the eigen- 
channel method and then the atomic multi-channel T- 
matrix and the dipole transition matrix elements are 
readily obtained from Eqs (|10|l and H15|l . We get the re- 
flectivity pLL'{ka) at the photo-electron energies of 
the different channels a, needed in Eq. ^ by interpola- 
tion in k. Finally, we invert the matrix t^^ — p (Eq. (S)) 
and obtain the XAS cross section from Eq. 

By virtue of the separation between environment and 
absorber through the partitioning technique, the present 
implementation of the multi-channel MS method is nu- 
merically only little heavier than the standard (single- 
channel) MS method. Indeed, in the present application, 
the atomic multi-channel calculation (second step above) 
was an order of magnitude faster than the reflectivity cal- 
culation by the standard MS technique. 
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FIG. 3: X-ray absorption spectra at the Ca-L2,3-edge in 
bulk Ca. Curve (e) is the experimental spectrum taken from 
Ref. [TUl . In the theoretical spectra, (g) was obtained from 
the ground state potential, and the others with a (partially) 
screened potential with the screening factor indicated as (a). 
In all cases, the one-electron spectrum (without V) is shown 
with a dashed line and the multi-channel calculation (includ- 
ing V^) with a full line. 



V. RESULTS FOR THE Ca La a-EDGE 

FigureOlshows the I/2,3-edge absorption of bulk Ca cal- 
culated in different approximations, along with the exper- 
imental spectrum (e) taken from Ref. [TJl • The numbers 
in parenthesis indicate the value of the screening param- 
eter a of the core-hole potential Vc in Eq. (fTB|l . The spec- 
tra labeled (g) have been obtained with the ground state 
potential (i.e. Vc = 0). The spectra in full (dashed) lines 
have been calculated with (without) the multipole part 
of the electron-hole interaction V. For easy comparison 
of the line shapes, all spectra are aligned at threshold 
and normalized with respect to the height of their main 
peak. Note that before normalization, the intensity of the 
spectra without V (dashed lines) was considerably bigger 
than the corresponding spectra with V (full lines). The 
relative renormalization factors between the two types of 
spectra, that have been used in Fig. 01 are: 1.8 (g), 2.3 
(0), 3.3 (0.1), and 3.8 (0.15). 

Probably the most striking feature of the spectra in 
Fig. is the effect of the multipole part of the electron- 
hole interaction V: in all cases, it leads to a big transfer 
of spectral weight from the L3 edge (lower energy peak) 
to the L2-edge. The branching ratio thus changes from 



2 : 1 without V to somewhat less than 1:1, which is in 
good agreement with experiment. This spectral weight 
transfer comes from the mixing between the 2pi/2 and the 
2p3/2-hole states (which correspond, in the one-electron 
approximation, to the L2 and L3 edges, respectively). 
It is a genuine atomic multiplet effect which was first 
explained by Zaanen et al. As can be seen from a 
"vertical" comparison in Fig. |2| the choice of the core- 
hole potential Vc has only a minor effect on the branch- 
ing ratio, but it changes the line shape of the two edges 
individually, as it can be expected from a single-particle 
quantity. Going from (g) to (0), or increasing the param- 
eter a has the effect of shifting the peak positions of the 
two edges to lower energy and of reducing their width. 
It moreover leads an overall shift of the whole spectrum 
to lower energy. This shift, which is roughly 1 eV for 
(g)^(O), (0)^(0.1), and (0.1)^(0.15), is, however, not 
apparent from Fig. 13 because we have aligned the spectra 
at threshold. When comparing the spectra including V 
with the experimental one, it is clear that (g) and (0) 
have much too broad peaks. Moreover, their peak posi- 
tions relative to threshold are at too high energy, espe- 
cially for (g). Good agreement for both peak width and 
positions is obtained for spectra (0.1) and (0.15). The 
only disagreement is that these two theoretical spectra 
show a weak fine-structure which was not observed ex- 
perimentally. A possible explanation for this discrepancy 
is the presence of further broadening mechanisms, other 
than coupling to the band, which is included here by 
the multiple scattering of the photo-electron. Himpsel 
et al. suggested that the broadening might be due to 
strong auto-ionization. The discrepancy could, however, 
also reveal limitations of the present screening model, 
which neglects charge fluctuations. Let us note that our 
spectrum (g) looks identical with the one obtained by 
Schwitalla et al. within time-dependent local density 
approximation. This shows that their method does not 
take account of the monopole part of the electron-hole 
interaction. 

For a contrast to the metallic bulk Ca, we have ap- 
plied the method also to two insulating Ca compounds: 
CaO and CaF2. The results are shown in Fig. 01 As 
in Fig. the spectra have been normalized and aligned 
at threshold. For the latter, the spectra (g) have been 
shifted by 3.5 eV to lower energy relative to the others in 
both compounds. The meaning of labels and line-styles 
is the same as in Fig. |21 The spectra (g), which corre- 
spond to a total neglect of the core-hole and the multipole 
terms V, are again completely at odds with the exper- 
imental spectrum (e). When using a screened core- hole 
potential Vc with a — 0.1, but still neglecting the mul- 
tipole terms V [dashed line (0.1)], the spectra consist of 
four narrow lines (the finite width comes entirely from 
the added Lorentzian-|-Gaussian broadening). The split- 
ting of the L3 and L2 peaks into two doublets is due to 
a strong ligand field effect. By symmetry resolved MS 
calculations, we have checked that for CaO, the lower 
(higher) energy peaks correspond to t2g (eg) symmetry 
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FIG. 4: X-ray absorption spectra at the Ca-L2,3-edge in CaO 
and CaF2 . Labels and line-styles have the same meaning as in 
Fig. El (s) experiment, (g) ground state potential, (0.1) Vc 
with Q=0.1. A full (dashed) line corresponds to a calculation 
with (without) V . 

states in the Oh point group. In CaF2 the order between 
t2g and Cg peaks is reversed. These spectra are, however, 
still very different from the experimental ones. When fi- 
nally also the multipole part of the interaction V is taken 
into account [full line (0.1)], very good agreement with 
experiment is obtained for both compounds. It should 
be noted that Himpsel et al. who used an atomic 
crystal field model, could also get very good agreement 
with experiment. However, in that work the crystal field 
is introduced empirically and its parameter values are 
adjusted to experiment. 

VI. CONCLUSIONS 

In summary, we have presented a method for X-ray ab- 
sorption in condensed matter where single-electron fea- 
tures are described in the MS approach, while local multi- 
electron effects are taken into account in a configuration 
interaction scheme. The novel features of the method are 
the multi-channel extension of MS theory and the use of 
an R-matrix technique in condensed matter. 

The method has been apphed to the Ca-L2,3 edge 
absorption of several Ca systems. The electron-hole 
Coulomb interaction was divided into its monopole and 
its (higher order) multipole part. The latter, which is 
responsible for the non-statistical L3:L2 branching ratio, 
was taken unscreened. We showed that no rescaling for 
this part is needed in our method in contrast to single 
configuration multiplet calculations. For the monopole 
term, a mixture between an unscreened and a statically 
screened core-hole potential was applied. A mixing fac- 
tor of about 10% yields line shapes in good agreement 
with experiment in all cases. 



Non-local correlation effects such as charge transfer ex- 
citations have been neglected in the present work. Let us 
mention, however, that such effects can, in principle, be 
included when the R-matrix reaction volume is extended 
from a single atom to a small cluster of atoms around the 
absorber. Compared to recent approaches based on time 
dependent density functional theory IT^ , we believe 
that the present, configuration interaction based method 
provides more insight in the correlation mechanisms at 
play. Moreover, the present approach can easily be ap- 
plied to problems where the applicability of TD-DFT has 
yet to be proved, namely core-level spectroscopies that 
involve more than one hole (such as Auger processes) 
or open 4/-shells. Compared to atomic multiplet meth- 
ods 10], the present approach does not rely on adjustable 
crystal field parameters. Instead, ligand field and band 
effects are described in an ah initio manner through MS 
theory. 
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APPENDIX 

Derivation of X-ray absorption cross section formula 
in multi-channel multiple scattering theory 

In this section we derive X-ray absorption cross sec- 
tion formula within the multi-channel multiple scattering 
method. Most of the results of this section are not new, 
but can be found as special cases of the more general 
derivation given in Ref. 0. Nevertheless, we think it is 
worth including this section for the convenience of the 
reader, because the derivation is simpler than the one 
in Ref. Q and it lends itself better to the form of the 
A'^-electron wave functions used in the present work. 

We start from the general multi-electron formula for 
the total optical absorption cross section in the dipole 
approximation P]: 

a = A^^au 5{Ef - E, - u) (17) 

/ 

Here, g and f are A^-electron wave functions, for ini- 
tial (=ground) and final state, respectively, in the ab- 
sorption process of a photon with energy ui. In case of 
degenerate ground states, a sum over g is understood. 
D = e • J2iLi is the dipole operator, a = the fine 
structure constant. We use the units: h = 1, Bohr radius 
for length, Rydberg for energy. Thus i^kin — k'^, — 2. 

For the ground state, we explicitly take into account 
only localized electrons of the absorbing atom. Thus 
we assume that the ground state wave function '^g is 
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confined to the atomic sphere of the absorber with ra- 
dius Tq: . . . xjv) = if 3i : Tj > rp. As for the 
final state wave function 5*/ , we assume that — 1 elec- 
trons remain in localized orbitals and at most one elec- 
tron (the "photo-electron" ) is promoted to a continuum 
orbital. Wc chose boundary conditions such that in the 
remote past, the photo-electron is free, i.e. its eigen- 
states are plane waves exp(ikx) times a spin function 
Xs(o') = Sscr- The rest system is in one of the eigen- 
states ^a{xi . . . xj^i^i) of the N—1 electron Hamiltonian 
with a core hole: H^^^^^ — Ea^a- Thus, the "incom- 
ing part" of is given by x exp(ikx)xs(cr). In multi- 
channel scattering theory, not only elastic, but also in- 
elastic scattering processes are taken into account, which 
correspond to excitations $q — *■ $/3. In the present ap- 
proach, these excitations are limited to atomic-like ones, 
such as multiplet excitations, due to the local character 
of <i>Q. (Note that this is in contrast to the more general 
theory in Ref. 1]). By expanding the scattered part of 
^' f over the eigenfunctions <i>Q , we can write 



^.("ks) ^ exp(zkx)x.((T) 



P 



fa (xct). 



Here /^"''*^(x(t) behaves asymptotically (r — * oo) like a 
purely outgoing spherical wave. In the above form of 
anti-symmetrization between the photo-electron and the 
A''— 1 other electrons has been disregarded. We indeed ne- 
glect anti-symmetrization for the "outside solution" , i.e. 
when the photo-electron is outside the atomic sphere of 
the absorber. For the solution inside the atomic sphere, 
however, anti-symmetrization between all electrons is 
correctly taken into account through the eigen-channel 
method (see main text). Note that in this work we have, 
for simplicity, assumed the "mufhn-tin" or more pre- 
cisely atomic sphere approximation for the one-electron 
potential i.e. the atomic cells are replaced by space fill- 
ing spheres with spherically symmetric potential inside. 
The difficulties of multiple-scattering theory arising from 
non-muffin-tin potentials are essentially independent of 
the electron correlation problem we are dealing with here. 
The present multi-channel approach could easily be gen- 
eralized to non-muffin-tin multiple scattering methods in 
which the muffin-tin spheres are replaced by space-filling 
atomic cells. The main change would consist in calculat- 
ing the R-matrix for a sphere surrounding the atomic cell 
and where the potential in the so-called "moon-region" 
(the space outside the cell and inside the sphere) has been 
put to zero. 

With the final state quantum numbers aks, the sum 
in Eq. 113 becomes Eaa/H^- We have / = 
/ dk de^/e, where e = fc^ is the kinetic energy 
of the photo-electron. This yields 



auj 



dkr, 



where fc^ — Eg + uj — Ea from energy conservation. 
It is convenient to work in an angular momentum ba- 



sis, i.e. to use spherical rather than plane waves. We 
have /*|k)(k| = IGtt^ |fci)(/fcL|, where (x|fcL) = 
ji{kr)YL{x) = JL{kx.). Here, ji are the usual spherical 
Bessel functions and Yj^ are spherical harmonics. The 
cross section now becomes 



(18) 



Here vl/("^'') is the scattering state that evolves from the 
incoming wave 



*'"^ = *a JL(fcaX)Xs(a) 



(19) 



Following standard multiple scattering theory, we write 
the scattered part of the wave as a sum of outgoing ir- 
regular waves from all the centers i 



In the following we consider points where the photo- 
electron coordinate x lies outside any muffin-tin sphere 
(or atomic cell). For such points the potential is zero and 
we have 



olLs 



(20) 



Here x^ = x — and HL{kx) = {kr)YL{x) where 
Ji^ = ji + ini is a Hankel and n/ a spherical Neumann 
function. As indicated by the superscript 0, the ampli- 
tudes -Bfais depend on the quantum numbers of ^i'""^, 
which we shall denote aoLoso from now on. We use the 
well-known re-expansion theorems: 

JL{kx,)=J2jL'ikx^)A'i,L{k) (21) 
L' 

-zi?L(fcx,) = E JL'{kxr)G%L{k) (22) 

L' 

where A^,^ and G^,^ are the real space KKR struc- 
ture constants. 01 Developing 5''"'^ and the 5'^'^'s around 
some given center i and using equations (|2l() and H22|l . 
respectively, yields: 

* = E {JL{kaXi}A",^Ls - ikc,HL{kaX^)B^^Ls} 

aLs 

(23) 

with 

^Lls — ^aaoSsso^LLoi^a) + ka E ^ LL' i^a) B^aL' s : 

(24) 

where the usual convention G^^, = has been used. 

Next we express the exciting wave amplitudes A*^^]^^ 
in terms of the scattered wave amplitudes B^^^^ at the 
same site i through the inverse atomic scattering matri- 

i^i )aLs,a' L' s' 8-8 

A-iaLs = E ^)aLs,a'L's'Bi^,i^,g, . (25) 
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This holds by definition of the ti-matrices, and rehes only 
on the most basic assumption of multiple scattering the- 
ory, namely that the potential can be written as a sum of 
atomic cell potentials. Note that there is no restriction 
on the form of the atomic potentials, which may, as it is 
the case for the absorber potential in the present work, 
include non-local and correlation effects. (Note, how- 
ever, that in the present approach the calculation of this 
complicated potential is avoided by virtue of the eigen- 
channel method.) Assuming the atomic t-matrices to be 
known, we may use Eq. (|25|l to eliminate the A'^^Ls'^ 
Eq. if^ . and then solve for the B'^^^^'s. This yields 

BiaLs — ^ '''aLs,aoi'so^i'Lo(^"o) (26) 

where r is the (multi-channel) scattering path operator 
which is defined by its matrix inverse: 

{t-^Y^^, = S,,{tr\r' - S^o.'ko.G'l^,{k^)Sss' , (27) 

where we have introduced the collective index T = aLs. 
Equations H23H27(I are the generalized multiple scattering 
equations. 

We shall proceed by calculating the X-ray absorption 
cross section from an atom placed at the origin = 0. 
Using Eq. ((SsJ, the wave function in Eq. around 
site i = (index suppressed) reads 

*ro -^3'r^rr'(r)AS^:? (28) 

where 

^t = <^o.Yl{x)xs{(t) (29) 

and 

ZrT'{r)/r = ]i{Kr){r^)rT' - iKhi{Kr)5rT' ■ (30) 

We recall that Eq. (|23|l or (|28|l is valid only in the space 
outside atomic spheres. For the region inside the atomic 



sphere of the absorber, we may write 

*ro-E*rSr° (31) 
r 

where ^j? is a solution of the Schrodinger equation in- 
side the atomic sphere, that matches smoothly onto the 
outside wave function $p'Zp/p (r)/r. 

Putting this into the absorption cross section formula, 
Eq. lfTH|) . we obtain 

TT'Ta 

(32) 

Note that the restriction to ^f'" in the calculation of the 
matrix elements is valid since we have assumed that g 
vanishes outside the atomic cell. A further simplification 
of the formula can be achieved if we use the optical theo- 
rem, whose validity in the multi-channel case was proved 
in Ref. JJ: 

Y.kv.B^rB]^r = -j^{r-r^y^^, (33) 

To 

If we moreover introduce the notation Mp = (5'p |_D|^g) 
we finally obtain: 



a = -Ai^auj X ^ I X] ^^r^rr' Mv | (34) 

In this form, the cross section formula reads exactly as 
the well-known one-particle expression (see e.g. Ref. 21]). 
The fundamental difference is that the quantum numbers 
r contain internal degrees of freedom of the absorbing 
atom (channels a), which in the present case correspond 
to different multi-electron states. 
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